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Abstract. In this paper we derive general expressions for few-electron spectral 
functions of the one-dimensional Hubbard model for values of the excitation energy 
in the vicinity of the M th upper-Hubbard band lower limit. Here M = 1,2, ... is the 
rotated-electron double occupation, which vanishes for the ground state and is a good 
quantum number for all values of the on-site Coulomb repulsion U . Our studies rely on 
a combination of symmetries of the model with a recent finite-energy holon and spinon 
description of the quantum problem. We apply our general scheme to the one-electron 
addition spectral function, dynamical structure factor, and spin singlet Cooper pair 
addition spectral function. Our results provide physically interesting information about 
the finite-energy spectral properties of the many-electron one-dimensional quantum 
liquid. 
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1. INTRODUCTION 

Recent experimental studies of quasi-one-dimensional (ID) materials observed unusual 
finite-energy/frequency spectral properties, which are far from being well understood 
[JJ El El HI IHj- Several of these experimental studies reveal the occurrence of 
charge-spin separation in terms of independent holon and spinon excitation modes 
[2 IU El HI El EH CD]- For values of the excitation energy larger than the transfer 
integrals associated with electron hopping between the chains, the ID Hubbard model 
[TT] is expected to provide a good description of the physics of these materials jH E]. 
Unfortunately, an accurate determination of the spectral properties for finite values of 
the excitation energy and of the Coulomb repulsion is until now still lacking. Most 
accurate results correspond to the limit of infinite on-site Coulomb repulsion ^2] where 
the Bethe-ansatz [TT] wave function is easier to handle, yet the solution of the problem 
remains complex in this limit. 

In this paper we combine the holon, spinon, and pseudoparticle accurate description 
recently introduced and studied in Refs. Ej and the pseudofermion representation 
very recently introduced in Ref. [Hj with symmetries of the ID Hubbard model to 
derive general expressions for few-electron spectral functions in the vicinity of the lower 
limit of the upper Hubbard bands. For finite values of the excitation energy and on-site 
repulsion U there are not many previous studies of few-electron spectral function weight 
distributions. The interesting studies of the one-electron spectral functions presented in 
Ref. ^2] refer to the limit of infinite U. There are numerical studies of these functions 
for U ~ At [JJJ]- Moreover, the weight distribution of the real part of the frequency 
dependent optical conductivity in the vicinity of the optical pseudogap was previously 
studied in Refs. jHHTJJj. The very recent results presented in Ref. [Ej refer to the one- 
electron spectral-function weight-distribution in the vicinity of singular branch lines. 
Some of these lines cross the upper-Hubbard band lower-limit points considered in 
this paper. However, our results are complementary of those presented in Ref. |lUj . 
which were obtained by use of the method of Ref. ^H] . Indeed, our weight-distribution 
expressions are valid when one approaches these points along all finite-weight directions 
except those corresponding to the above-mention singular branch lines. The values of 
the critical exponents which control the weight distribution are different for these lines 
and for all the remaining directions of the (k, u;)-plane considered here. 

We evaluate expressions for the critical exponents that control the weight 
distribution of few-electron spectral functions for excitation energies corresponding to 
the vicinity of the M th upper Hubbard band lower-limit points. (The method used in 
this paper is different from that of Ref. [IB] , which was used in Ref. [E] in the study 
of the weight distribution in the vicinity of the singular branch lines.) Our studies rely 
on a symmetry which is specific for the model when defined in the reduced Hilbert 
subspace associated with the (k, c<j)-plane region in the vicinity of the above lower-limit 
points. Few-electron spectral functions can include several upper Hubbard bands, the 
M th band being spanned by excited states of rotated electron double occupation M > 0. 
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The lower Hubbard band corresponds to the spectral-weight distribution generated by 
excited states of rotated-electron double occupation M — 0. For one-electron and two- 
electron spectral functions only the first few upper bands have a significant amount of 
spectral weight [16J. Although our general expressions refer to the lower-limit points of 
all upper Hubbard bands, we apply our method only to the study of weight distributions 
corresponding to upper bands with a significant amount of spectral weight. For the one- 
electron addition spectral function and dynamical structure factor (and the spin singlet 
Cooper pair addition spectral function) only the M = 1 first Hubbard band has (and 
both the M = 1 first and M = 2 second upper Hubbard bands have) a significant 
amount of weight. 

Our motivation in studying the finite-energy one-electron spectral function is that 
our method provides the exponents of the singular spectral features beyond the branch 
lines considered in Ref. jTOj. The latter lines are observed in quasi- ID materials |H| ll()j. 
Thus, the study of other finite-energy singular spectral features is of interest for the 
further understanding of the unusual spectral properties observed in these materials. 
Both for the dynamical structure factor and other spectral functions we do not consider 
the low-energy features because those can be investigated by standard two-component 
conformal-field theory. Our aim in studying the Cooper-pair spectral function is to 
clarify whether there are singular spectral features at the upper-Hubbard bands lower 
limit. Indeed, in the limit of small on-site repulsion and for electronic densities very 
close to one, such singular features would appear at low energy and could lead to a 
superconductivity instability in a system of weakly coupled Hubbard chains, as further 
discussed in later sections. 

The paper is organized as follows: In Sec. II we summarize the ID Hubbard model. 
In Sec. Ill we introduce the holon-spinon Hamiltonian and study the ID Hubbard model 
spectrum in the vicinity of the lower limit of the upper Hubbard bands. In Sec. IV 
we use the holon and spinon conservation laws and other symmetries of the model in 
the evaluation of general expressions for few-electron spectral functions. Moreover, we 
apply our method to the study of the one-electron addition spectral function, dynamical 
structure factor, and singlet Cooper pair addition spectral function. Finally, Sec. V 
contains the concluding remarks. 

2. THE ID HUBBARD MODEL 

In a chemical potential fi and magnetic field H the ID Hubbard Hamiltonian can be 
written as, 

H = H S o(4) + ^2 ^ a > 

a=c, s 

N a N a 

H SO (a) =-*E E [cl ff c j+ i, a + 4+i,* c j,<t\ + U Et%,T - 1/2] [nj,i - 1/2] 

J'=l <T=U j = l 
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where c] a and Cj )Cr are the spin a electron creation and annihilation operators at site 
j, respectively. The operator h^ a = c] )0 .c J)(T counts the number of spin a electrons at 
real-space lattice site j = 1, N a . The number of lattice sites N a is even and large and 
N a /2 is odd. We consider periodic boundary conditions. Moreover, \i c = 2/i, // s = 2ft H, 
Ho is the Bohr magneton, and S c z = -(l/2)[N a - N] and S s z = -(1/2) [iV T - N L ] are 
the diagonal generators of the //-spin and spin SU (2) algebras |17J , respectively. The 
operators N = N a and N a = J2j a count the number of electrons and spin a 
electrons, respectively. The Hamiltonian HsoU) defined in Eq. (0) commutes with the 
six generators of the 77-spin and spin algebras ^7]. The Bethe-ansatz solvability of the 
ID Hubbard model is restricted to the Hilbert subspace spanned by the lowest-weight 
states (LWSs) of the 77 spin and spin algebras, i.e. such that S a = —S° jH] for a = c 
and a = s, respectively. Here S c (and S s ) denotes the ?7-spin (and spin) value. The 
momentum operator reads, 

p= E (2) 

where the spin a momentum distribution operator is given by N a (k) = c\ a Ck, a and 
the operator c\ a (and Cfc )CT ) creates (and annihilates) a spin a electron at momentum k. 
The momentum operator (j2J) commutes with the Hamiltonians of Eq. $1} . 

There are spin-up electrons and N± spin-down electrons in the chain of iV a sites, 
lattice constant a, and length L = [N a a] associated with the model (JTJ). Throughout 
this paper we use units of Planck constant one and of lattice spacing a = 1 and 
denote the electronic charge by — e. The Fermi momenta read kp a = 7m a and 
kp = ±[fci?| + kp\\/2 = irn/2, where n a = N a /N a and n = N/N a . The electronic density 
can be written as n = nj + and the spin density is given by m = n-j- — n^. In general 
we consider electronic densities n and spin densities m in the domains < n < 1 and 
< m < n, respectively. Our general correlation-function expressions are derived for the 
metallic phase of electronic densities < n < 1 and spin densities < m < n. However, 
taking the limit n — > 1 in some of these expressions leads to correct expressions for the 
n — 1 Mott-Hubbard insulator phase. Also the m = spectral function expressions 
can in general be obtained by taking the limit m — > in our general expressions. In 
our applications to the study of specific spectral functions we consider such m = 
expressions only. 

The concept of rotated electron [HI 13 IE] is associated with a unitary transformation 
introduced in Ref . [TH] . For such rotated electrons, double occupation is a good quantum 
number for all values of the on-site Coulombian repulsion U. The electrons that occur 
in the ID Hubbard model (JTJ) are defined by c'j a , while the rotated electron operator 
cj CT is given by, cj j(7 = V~^(U/t) Cj a V(U/t), where V(U/t) is the electron - rotated- 
electron unitary operator defined by Eqs. (10)-(12) of Ref. |H]. As a result of the 
Hilbert-space electron - rotated-electron unitary rotation, all energy eigenstates of the 
model are described in terms of occupancy configurations of s c = 1/2 and a c = ±1/2 
holons, s s = 1/2 and a s = ±1/2 spinons, and c pseudoparticles Here we denoted the 
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77-spin and spin projections of the quantum objects by a c and a s , respectively, whereas 
their //-spin and spin values are denoted by s c and s s , respectively. In this paper we 
call the holons and spinons according to their values of a c = ±1/2 and a s = ±1/2, 
respectively. We denote by M^ ±lj / 2 the number of ±1/2 holons (a = c) or ±1/2 
spinons (a = s). The value of the number of rotated-electron doubly occupied and 
unoccupied sites (and spin-down and spin-up rotated-electron singly occupied sites) 
equals that of the holon numbers M C) _i/2 and M c +1 / 2 (and spinon numbers M S! _i/ 2 
and M S)+ i/ 2 ), respectively. Thus, rotated-electron double occupation M is such that 
M = M Cj _ 1/2 . Moreover, M a = M Qj+1/2 + Af a> _i/ 2 and N c = M s , where N c is the 
number of c pseudoparticles. The ±1/2 holons and c pseudoparticles carry charge ±2e 
and — e, respectively, whereas the spinons have no charge degrees of freedom. The cv 
pseudoparticles (and sv pseudoparticles) associated with Takahasi's charge (and spin) 
ideal string excitations of length v |11[IH] are 7/-spin singlet 2v- holon (and spin singlet 2v- 
spinon) composite quantum objects |H10IE]- The ±1/2 holons (and ±1/2 spinons) which 
are not part of such 2i/-holon composite cv pseudoparticles (and 2i/-spinon composite 
sv pseudoparticles) are called ±1/2 Yang holons (and ±1/2 HL spinons) 0. In the 
designations HL spinon and Yang holon, HL stands for Heilmann and Lieb and Yang 
refers to C. N. Yang, respectively, who are the authors of the papers of Ref. |17j . 
We denote by N au the number of composite av pseudoparticles belonging to branches 
a = c, s and v — 1,2, .... We call L a ,±i/ 2 the number of ±1/2 Yang holons (a = c) or 
±1/2 HL spinons (ct — s). Note that M a ,±i/2 — -^a, ±1/2 ± S^Li v N au . 

Below we are mostly interested in excited states without HL spinons and sv 
pseudoparticles such that v > 1. Therefore, for simplicity often we replace the 
pseudoparticle branch index si by s. Thus, throughout this paper the notations si 
and s are equivalent. For the subspaces spanned by the ground state and excited states 
generated from it by a finite number of pseudoparticle and/or Yang holon processes, 
the c, s, and cl pseudoparticles carry band momentum q such that \q\ < g° = n, 
I <?| < <7s = kp, and \ q\ < q^ = [n — 2k F ], respectively. This refers to electronic densities 
< n < 1 and spin density m = 0. For these densities the ground state is such 
that there are no —1/2 Yang holons, the cl and s pseudoparticle bands are empty and 
filled, respectively, and the c pseudoparticles occupy < | q\ < q Fc — 2kp and thus 
leave 2k F < \ q\ < 7r empty |E1IZ|- Following the notation of Ref. [Ij, we call an electron 
ensemble space a Hilbert subspace spanned by all states with fixed values for the iVj and 
iVj electron numbers. Furthermore, we call a CPHS ensemble space, where CPHS stands 
for c pseudoparticle, —1/2 holon, and —1/2 spinon, a Hilbert subspace spanned by all 
states with fixed values for the numbers iV c , M Cj _i/ 2 , and M Sj _i/ 2 . A CPHS ensemble 
subspace is a Hilbert subspace spanned by all states with fixed values for the numbers 
N c , L c _!/ 2 , and L Sj _i/ 2 and for the sets of numbers {N cv } and {N su } corresponding to 
v = 1, 2, ... branches. We note that one does not need to provide the values of M c +1 / 2 
and M Sj+ i/2 in order to specify a CPHS ensemble space, since these numbers are not 
independent |7j. Also the numbers L Cj+ i/ 2 and L Sj+ i/ 2 are not independent and one 
does not need to provide these values in order to specify a CPHS ensemble subspace. 
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Finally, let us introduce the useful quantum number i = sgn(q)l = ±1 which refers 
to the number of right pseudoparticle movers (l = +1) and left pseudoparticle movers 
(a = —1). The numbers N C>L of c pseudoparticles and N aVtL of av pseudoparticles of i 
character are good quantum numbers. We thus introduce the c pseudoparticle and av 
pseudoparticle current numbers, 

■Je^EWJVc,*; J» = 5EW^- ( 3 ) 
z t=±i z t=±i 

The numbers N C)L and iV^t can be expressed as N C)i = N c /2 + lJ c and N av>l = 
N au /2 + iJ av . 

Each CPHS ensemble Hilbert subspace characterized by fixed values for the sets 
of numbers N c , {N au } such that a = c, s and v = 1,2,3,..., and {L a _i/ 2 } such that 
a = c, s contains different subspaces with different values for the sets of current numbers 
J c and { J au } such that a = c, s and v — 1, 2, .... According to the notation of Ref. [S], 
we call these subspaces J- CPHS ensemble subspaces. 



3. THE HOLON-SPINON HAMILTONIAN AND THE UPPER 
HUBBARD BANDS LOWER LIMIT SPECTRUM 

The correlation function expressions are determined by transitions from the ground state 
to excited states. Let us consider that the initial ground state belongs to a canonical 
ensemble space associated with densities such that < n < 1 and < m < n. Often 
we consider such ground states. In the application of our results to the evaluation of 
specific spectral function expressions we are mostly interested in m = initial ground 
states. Fortunately, we are able to derive the corresponding m = expressions from the 
general correlation functions expressions obtained for < m < n. Let us introduce the 
Hamiltonians Hql and Hhs such that, 

Hgl = H — H HS ; 

oo 

H HS = 2fi M c , _ 1/2 + 2/i H [M s , _ 1/2 - N sl ] + £ e° sl/ (0) N 8V . (4) 

Here H is the ID Hubbard model Hhs is the holon-spinon Hamiltonian, and in 
H GL the letters GL stand for gap less. Indeed, the energy spectrum of the Hamiltonian 
Hgl of Eq. (jlj is gapless. Moreover, M a ,-i/2 is the —1/2 holon (a = c) and —1/2 
spinon (a = s) number operator given in Eqs. (24) and (25) of Ref. jH], N S u is the sv 
pseudoparticle number operator, and e° sv (q) is the sv pseudoparticle energy band defined 
in Refs. [01 13 El- For m — > and v > 1 this energy band is such that e° sv (q) = 0. 

Since the numbers of —1/2 holons, —1/2 spinons, and sv pseudoparticles are good 
quantum numbers, the three Hamiltonians of Eq. (jlj) commute with each other. Thus, 
these Hamiltonians have the same energy eigenstates, only the corresponding energy 
eigenvalues being in general different. Let us consider an energy- E eigenstate of the 
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ID Hubbard model. Such a state is also an energy eigenstate of the holon-spinon 
Hamiltonian H HS whose energy eigenvalue ujhs reads, 

oo 

lu hs = 2fM M c> _ 1/2 + 2/i H [M Sj _ 1/2 - N sl ] + £ 4(0) N w . (5) 

It is also an energy eigenstate of the Hamiltonian Hql of Eq. (J1J of eigenvalue Eql = 
[E — Ujis]- Analysis of the finite-size corrections of the energy AE GL = [Eq L — Eqs] ~ 
of all low-energy states of the Hamiltonian Hql reveals that the ground state of the ID 
Hubbard model is also the state of minimal energy AE GL = [E GL — E GS ] = for such a 
Hamiltonian. Indeed, for all canonical ensemble spaces the latter Hamiltonian and the 
ID Hubbard model (JIJ have the same ground state. The general spectrum of interest 
for the problem of the few-electron spectral functions of the ID Hubbard model for 
excitation energies uj = [E—Eqs] = uhs + AE G l such that AE G l = [E G l — E G s is small 
can be written as the sum of the finite energy ojhs plus the gapless contribution AE GL 
expressed in terms of the pseudoparticle energy bands [HE]- This gapless contribution 
is the excitation energy of the same states relative to the Hamiltonian H G l of Eq. 

The ground-state normal-ordered expression of an operator O is defined as : O : = 
O — (GS\ 0\ GS), where | GS) is the ground state. Let \ex) be any excited state. The 
corresponding deviation AO is defined as AO = (ex\ \ ex) — (GS\ 0\ GS). From now 
on we will deal mostly with the two Hamiltonians and H G l of Eq. (0J) while acting 
in the reduced Hilbert subspaces spanned by excited states with small excitation energy 
AE G l for the Hamiltonian H G l- Thus, it is useful to introduce the ground-state normal- 
ordered Hamiltonians : Hgl '■ an d : H : such that, 

:H GL :=:H:-H HS , (6) 

where H is the ID Hubbard model (JTJ). (Note that since the energy eigenvalue u>hs © 
of the holon-spinon Hamiltonian Hhs vanishes for the ground state, its ground-state 
normal-ordered expression is such that : Hhs '■= Hhs-) While the Hamiltonian : H G i : 
of Eq. (jSJ describes the gapless part of the excitations, the holon-spinon Hamiltonian 
Hhs of Eqs. (@J) and (JHJ) controls the finite-energy physics. Creation of a —1/2 ho Ion 
requires a minimal amount of excitation energy 2/x. Creation of a —1/2 spinon (except 
those which are part of si pseudoparticles) requires a minimal amount of excitation 
energy 2u,$H . Creation of a sv pseudoparticle of bare-momentum q and belonging to a 
v > 1 branch involves creation of v —1/2 spinons and is associated with an energy given 
by 2u fj, H+e^ u (q). Creation of a cv pseudoparticle of bare-momentum q involves creation 
of v —1/2 holons and requires an amount of energy 2ujj, + e® u (q). (The pseudoparticle 
energy bands e° sv (q) and e° cv (q) are studied in Refs. jHllZl-) 

Although our expressions refer to general values of spin density m such that 
< m < n, in the case of applications to few-electron spectral functions we are most 
interested in zero-magnetization. Thus, for simplicity we consider excited states without 
— 1/2 HL spinons and without sv pseudoparticles belonging to v > 1 branches. We 
emphasize that finite-energy states with finite occupancies for these quantum objects 
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play an important role only for finite values of the spin density. One arrives to the same 
final m = expressions from our general < m < n scheme regardless we consider 
excited states with or without —1/2 HL spinons and sv pseudoparticles belonging to 
v > 1 branches. 

For m —>■ intial ground states only excited states with finite —1/2 ho Ion 
occupancies have a gapped energy spectrum. These states have a minimal finite 
excitation energy given by, 

oo 

uj hs = E u M Cj _i /2 = E u [L Ci _i /2 + Yl v N cu\ > ( 7 ) 

Here the energy E u = 2fi is defined by Eq. (107) of Ref. |S|. In the limit m —>■ 
and for < n < 1 it is an increasing function of the on-site repulsion U such that 
E u = 4t cos(7rn/2) for U/t — > and E u = U + 4t cos(7rn) for U/t — > oo. For any value of 
U/t it is a decreasing function of the electronic density n such that E u = U + 4t as n — > 
and ^ approaches the value of the Mott-Hubbard gap Emh as n — > 1. For values of 
rotated-electron double occupation M such that M > 0, the energy luhs = M E u of Eq. 
(JZJ) is finite and corresponds to the lower limit of the M th upper Hubbard band. 

Transitions from a ground state to excited states belonging to a given J-CPHS 
ensemble subspace can be labelled by the set of deviation numbers and deviation 
current numbers {AN a }, {AJ a }, {AN CU } = {N cu }, {AJ CU } = {J cu }, and {AL Ci _ 1/2 } = 
{Lc-1/2} where a = c, s and v — 1, 2,.... Since there are no cv pseudoparticles and 
—1/2 Yang holons in the ground state, below we replace the deviations AN CU , AJ CU , 
and AL C) _!/ 2 by the corresponding number and current number values N cu , J cu , and 
L c ,-i/2, respectively. All states belonging to the same J-CPHS ensemble subspace have 
the same number of —1/2 holons and thus the same energy eigenvalue ujhs — M E u 
relative to the holon-spinon Hamiltonian Hrs of Eq. (£Q). We are mostly interested in 
the reduced J-CPHS ensemble subspaces. Such subspaces are the part of the general 
J-CPHS ensemble subspaces which is spanned by states of excitation energy u (relative 
to the ID Hubbard model) such that AEgl = {uj — M E u ) is small. Transitions from a 
ground state to these states correspond to transitions to the lower limit of the M th upper 
Hubbard band. Importantly, the cv pseudoparticle creation processes which generate 
the M Ci _!/2 = M = 1,2, ... excited states of vanishing energy AEgl ~ and M > 
excitation energy AE = u ~ M E u for the ID Hubbard model involve cv pseudoparticles 
of bare-momentum value q ~ ±5^, = ±[ir — 2kp] only. In addition, these states can 
involve —1/2 Yang holons creation and low-energy c and s pseudoparticle creation 
and/or annihilation and particle-hole pseudoparticle processes in the vicinity of the 
corresponding c and s pseudoparticle Fermi points. Moreover, since we are considering 
final states with no —1/2 HL spinons and no sv pseudoparticles belonging to v > 1 
branches, the number deviations and numbers AiV c , AN S , L c _i/ 2 , and {N cv } where 
v — 1, 2,... obey the sum rules (51) and (52) of Ref. [7] with AL St _i/ 2 = and 
AN sl/ = for v > 1 (In that reference N sl/ is denoted by N SjU .) 

When defined in a reduced J-CPHS ensemble subspace, the momentum operator 
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(J2J) can be written in normal order relative to the ground state as follows, 
: P := P + : P GL : ; 

oo 

Po = E 9f*2 : J a :+J2 2 k F 2J cu + 7iM C: _ 1/2 . (8) 



i/=i 



Here Pq is the momentum associated with the holon-spinon Hamiltonian Hhs of Eq. 
(UJ). Moreover, the ground-state normal-ordered operator : P G l : of Eq. (JSJ) is the 
momentum operator associated with the Hamiltonian : Hql '■ of Eq. (JSJ) and reads, 



: Pgl '■— -rr E : ^ a : : ^ a : — ^ a > c E ^< 



CI' 



p, 



ph j 



jv a q=c,si=±1 

The operator P p h is such that A^ h t counts the number N%r L of momentum t[27r/A^ a ] 
elementary particle-hole pseudoparticle processes around the a = c and a = s Fermi 
points iq Fc = i2kp and iq Fs = ikpi, respectively. 

The momentum operators involved in Eq. (JSJ) commute with each other and 
with the three Hamiltonians appearing in Eq. (JHJ). Thus, the energy eigenstates are 
also eigenstates of the latter momentum operators. The corresponding momentum 
eigenvalues are such that, 

oo 

AP = k l M +AP GL ; k l M = £ g° Q 2AJ a +]T2A; F 2J cl/ +7rM Ci _ 1/2 .(iO) 

a=c, s u=l 

The index I labels the different momentum values k l M occurring for the same value of 
rotated-electron double occupation M, as illustrated in the applications provided below. 
The momentum APql of Eq. (TTUJ) reads, 

2tt 



AN C \AJ C - ^ J CIJ ] + AN S AJ S + ^2 Y. 1 N * 

u=l L =±la=c,s 
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The current number deviations AJ C , AJ S , and J cl/ = AJ CU of Eqs. (JTHJ) and (fTTjl are 
given by, 

A J c = ^ J modi ; AJ S = \ — — J modi ; 

J*, = \ E iA W- ( 12 ) 

The momentum values k l M defined in Eq. fllOj) play the same role for the J-CPHS 
ensemble subspaces as the energies u>hs — M E u for the CPHS ensemble spaces. While 
the straight horizontal line (k, uj = M E u ) defines the lower limit of the M th upper 
Hubbard band, the points (k = k l M , u = M E u ) of the same line define the edges of the 
M th upper Hubbard band, as further discussed in Sec. IV. We call 7i re d the reduced 
J-CPHS ensemble subspaces spanned by states of momentum k and excitation energy uj 
such that (k — k l M ) and {uj — M E u ) are small. These subspaces are characterized both 
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by specific values for the set of deviation numbers and numbers {AN a }, {N cu }, {A J a }, 
{J CI/ }, and {L Cj _!/ 2 } where a = c, s and v — 1, 2, ... and by small values for (k — k l M ) 
and (a; — M E u ). 

It is useful for our studies to find the specific form of the general energy spectrum 
AEql of the Hamiltonian : Hql '■ of Eq. (jOJ for the excited states which span a 
reduced subspace TL re d- That energy spectrum can be expressed in terms of a phase- 
shift momentum functional associated with pseudofermions introduced in Ref. 
According to the results of that reference, there is a canonical transformation that maps 
a pseudoparticles and cv pseudoparticles onto a pseudofermions and cv pseudofermions, 
respectively. While the a pseudoparticles and cv pseudoparticles carry bare-momentum 
q, the a pseudofermions carry momentum q = q + Q a (q)/N a and the cv pseudofermions 
carry momentum q = q + Q C v(q) /N a . The functionals Q a (q) and Q C v(q) are defined in 
Eq. (73) of Ref. |Hj. After some algebra one finds that for the above states AEql has 
the following leading-order 1/N a finite-size energy terms, 



2tt 



AE ol = W E E K + Kt +0(— ). (13) 



1 



a a=c,s i=±l ±y a 



Here v a stands for the light group velocity v a = v a {q° Fa ) where v a (q) = de a (q)/dq for 
a = c and a — s. Moreover, N% \ = 0, 1, 2, ... is the number of elementary particle- hole 
pseudoparticle processes around the i — ±1 Fermi points of the c and s pseudoparticle 
bands and 2 stands for the following functional, 



2 At. 



7^{ A <lFa,, + ^ — ^ — JJ 5 a = c,s; l = ±1, (14) 

where Q a {q)/N a is the above a pseudofermion momentum functional. Note that 
since q = q + Q a (q)/N a , the functional (jTJ| is proportional to the square of the 
a pseudofermion Fermi momentum deviation A(f Fa L = Aq Fa L + Q a (qFa J/-^o- m 
Appendix A it is found that for the excited states that span the reduced J-CPHS 
ensemble subspaces Ti re d the functional (|T4*j) can be written as, 

2K=[c e ^-2 £L +eic[AJc-E^]+C AJ -] • ( 15 ) 

a'=c,s v=\ 

Here the parameters ^ a , are defined by Eq. (|A.2|) of Appendix A. 

A property with a deep physical meaning is that the finite-size energy spectrum of 
the cv pseudoparticles vanishes as the limit q — > ±g^ = ±[n — 2k F ] is approached. In 
that limit the cv pseudoparticles become localized and non-interacting objects, as found 
in Ref. [7] and discussed in Appendix A. In contrast, creation of cv pseudoparticles at 
other bare-momentum values leads to finite-size cv energy corrections and to values of 
excitation energy u such that AEql — (to — M E u ) is finite. In that case the finite- 
energy problem cannot be mapped onto a low-energy conformal-field theory and our 
method does not apply. Behind this vanishing of the finite-size cv energy spectrum 
there is a symmetry which imposes that the v — 1/2 Yang ho Ion and cv pseudoparticle 
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energy and momentum spectra become the same as q — *> ±[ir — 2kp] for one of the 
bare-momentum values — [tt — 2kp] or +[tc — 2kp\. Creation of v —1/2 Yang holons 
contributes to the ID Hubbard model energy spectrum u> = ujhs + AEql through the 
finite-energy term u HS = M E u of Eq. (J7J) by an energy amount v E u , but does not 
lead to any contribution to the gapless spectrum AEql flEU) and thus to the value of 
the functional (fT3j) . Therefore, an excitation involving c and s pseudoparticle processes 
and creation of v —1/2 Yang holons leads to an energy spectrum of the form (fT^j) . 
whose functional 2 is given by (fT5J) with A J c — J2^Li Jcu replaced by A J c . Except 
for this difference, the cv pseudoparticle and v — 1/2 holon contributions to the energy 
spectrum u = ujhs + AEql and functional (fTKj) are the same. However, the role of the c 
pseudoparticle current number deviation shift — Y^=i Jcu is to introduce a counter term 
whose presence assures that two M > finite-energy excited states with the same c and s 
pseudoparticle occupancies and either one cv pseudoparticle such that u = M ot M — 1/2 
Yang holons, respectively, have the same momentum and energy spectrum. The same 
holds for general M > excited states such that M = YJ^° v N cu +L Ct _y 2 . For simplicity, 
let us consider that M — 1. In this case the current — Y^T=\ Jcu reads — J c i = ^1/2 and 
zero for creation of a cl pseudoparticle and a —1/2 Yang holon, respectively. According 
to Eq. ()12j) . for states with the same c and s pseudoparticle number deviations the 
value of the c current number deviation AJ C for the one-cl pseudoparticle states differs 
from the corresponding number of the one — 1/2 Yang holon states by ±1/2. Thus, 
for one of the two possible values of opposite sign of the current number — J cl = =Fl/2 
such difference is precisely cancelled. For that one-cl pseudoparticle state the current 
number AJ C — J cl equals the current number AJ C of the one — 1/2 Yang holon state. 

4. SPECTRAL FUNCTION EXPRESSIONS IN THE VICINITY OF THE 
UPPER-HUBBARD BANDS LOWER LIMITS 

In this section we derive general finite-energy expressions for few-electron spectral 
functions for excitation energy and momentum values in the vicinity of upper-Hubbard 
bands lower limits. Moreover, we apply these general expressions to the study of specific 
one-electron and two-electron functions. 

4.1. THE GENERAL FINITE-ENERGY SPECTRAL FUNCTION EXPRESSIONS 

Although for M E u > the energy spectrum u = M E u + AEql of the the ID 
Hubbard-model in the reduced subspace Tired refers to a finite-energy problem, the 
functionals 2A^, and 2AJ, appearing in Eq. (|T3j) equal the c and s primary-field 
dimensions of a low-energy two-component conformal-field theory [TTfl |2*U] . We note 
that the momentum spectrum (jllj) can be expressed in terms of these functionals as 
APgl = jf- H a =c,s S t =±i L \^a + -^afj • Furthermore, manipulation on the integral 
equations that define the phase shifts of the ^ a , and £° a , expressions given in Eq. 
(|A.2J) of Appendix B, reveals that those are entries of known matrices. The two 
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corresponding matrices are the transpose of the dressed charge matrix and the inverse of 
the transpose of the dressed charge matrix, respectively, of the low-energy ID Hubbard 
model conformal field theory ^01 ED]- This property allows the evaluation of the 
ID Hubbard-model few-electron spectral function expressions for momentum k and 
excitation energy to such that both (k — k l M ) and (u — M E u ) are small. For M E u > 
such a ID Hubbard model finite-energy problem cannot be solved by standard low- 
energy conformal field theory. 

Let (j)#(x, t) (and <j)^ L (x, t)) represent an one-electron or two-electron physical 
field and d = lp, p, ss refer for example to one-electron, charge, and singlet 
superconductivity, respectively. The time evolution of (j>#(x, t) (and <j)$ L {x, £)) is 
described by the ID Hubbard Hamiltonian : H : (and Hamiltonian : Hql '■ of Eq. 
(JEJ)). The space translations are described by the momentum operator : P : of Eq. (JHJ) 
(and : Pql '■ of Eq. (jHJ)). The asymptotic expression for the low-energy correlation 
function of the physical field (p^ L (x, t) can be obtained by use of conformal- field theory. 
Such an expression, combined with the relation between the three Hamiltonians of Eq. 
(JHJ and their symmetries, is used in Appendix B to show that the leading term in the 
asymptotic expansion of the corresponding finite-energy correlation function of the ID 
Hubbard model is of the following form, 

e -i[k l M x-M E u t] 

X*(x, t) = (GS\ Mx, t) <M0, 0) I GS) oc [I II TZ T^r ■ (16) 

a=c,s 4 =±1 \ X LV a^) a 

Here is the functional defined in Eq. (|15j) and v a stands for the light group velocity 
v a = v a (qp a ). Comparison with the low-energy correlation-function expression (jB.2|) of 
Appendix A for <p^ L (x, t) reveals that there is in expression (jl6|) an extra phase factor 
e -i[k M x-ME u t]_ j n S pi^ e f that, for values of the momentum k and excitation energy uj 
such that both [k — k l M ) and {uj — M E u ) are small the asymptotic of the finite-energy 
correlation function is of algebraic type. The non-interacting character of the —1/2 
Yang holons and cl pseudoparticles as the limit q — > ±gj?^ = ±[n — 2kp] is approached 
justifies that the creation of such objects is a finite-energy process that leads to the 
phase-factor e - l i kl M x - M E u t\ Q f eX p ress i on (jjgj) only. 

In order to derive the expressions of the few-electron spectral functions one needs 
the imaginary part of the correlation functions in the k and u plane. Let x^i^, 
be the Fourier transform of the function x${ x ) t) given in Eq. ((16)1 . The x^i^, 
expression is controlled by the same universal exponent when one approaches the point 
(k = k M , uj = M E u ) from all directions (k — > k l M , uj — > M E u ) in the finite spectral- 
weight region corresponding to the reduced subspace TC re d, except for the four lines such 
that {uj — M E u ) « ±v a (k — k l M ). (Two lines for each a = c, s pseudoparticle branch 
- in some cases only two of these four lines are inside the finite spectral- weight region.) 
These lines are associated with the slope at the point (k = k l M , uj = M E u ) of the 
general a branch lines studied for few-electron spectral functions in Refs. ^Hl • in 
the vicinity of these lines the uj dependence of the weight distribution is also in general 
of power-law type but is controlled by exponents different from that given below. Thus, 
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our studies are complementary to those of Refs. fOJEIj. For simplicity, we consider 
that the critical point is approached through the line defined by k = k l M and uj —>■ M E u . 
In this case we find the following expressions for x${k-> u), 

lm X 4k M , cx (u - ME u y* , (17) 

and 



ReX#(k l M , uj) oc (uj - 

oc -ln(u -ME U ); C* = 0, (18) 

where the exponent reads, 

o = -2+ E E 2 K, (19) 

a=c, s t=±l 

and the functional 2A L a is given in Eq. (|T3j) . Expression ()17|) corresponds to the leading 
term of an expansion in the small energy (uj — M E u ). The exponent is a rapidly 
increasing function of the number of excited pseudoparticles. The non-interacting 
character of the —1/2 Yang holons and q = cv pseudoparticles justifies that for 
m — > the U/t and n dependence of the exponent (|T^|) occurs through the parameter 
K p of Ref. [22] only, as confirmed by the exponent expressions found below. (This is 
not in general true for the branch-line exponents of Refs. |TTH ITH] , whose expressions 
involve the momentum-dependent phase shifts defined in Ref. jH]-) Importantly, when 
for one-electron (and two-electron) correlation functions the exponent tends to — 1 
(and —2) as U/t — > or U/t — > oo, expression (|T7j) is not valid. In this limit the 
correlation-function expression is such that, 

X#{k l M , w) oc — ; Ina.x#{±k l M , uj) oc 6(u - M E u ) , (20) 

where /3 is real and infinitesimal and thus the spectral function lmxd{±k l M , uj) is 8- 
f unction like. 

For simplicity, the general expression (j!7)l refers to the vertical line corresponding 
to to k = k l M and small values of (a; — M E u ), yet the exponent (JTHJ) controls the weight 
distribution for all other lines crossing the point (k = k l M , uj = M E u ) except for the a 
branch lines. Thus, although in our applications to the one-electron addition spectral 
function, dynamical structure factor, and singlet Cooper pair addition spectral function 
we use the k = k l M expression (fTTj) . we emphasize that similar expressions apply to the 
weight distribution associated with momentum values k and excitation energy uj such 
that {k — k l M ) and (lu — M E u ) are small. 

Except in particular limits of parameter space, our present method does not provide 
the values of the constants that multiply the weight-distribution power-law expressions. 
This justifies the use of the proportionally symbol oc in the general expression (fT7|) . 
Fortunately, there is a relation between the value of the critical exponents that control 
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the weight distributions in the vicinity of the different upper Hubbard band lower-limit 
points and the relative value of the corresponding power-law multiplicative constants. In 
general, the smallest is the critical exponent, the largest is the value of the corresponding 
multiplicative constant. Such a relation is confirmed by analysis of the results obtained 
for U/t — > oo by the method of Ref. |12j . The finite-[//t numerical results of Ref. ^3] 
and the small U / 1 results of Ref. also confirm that relation. Moreover, often the form 
of the weight distributions at U/t = also provides an useful boundary condition. Thus, 
our considerations for finite values of U/t concerning the relative amount of spectral 
weight located in the vicinity of the different upper Hubbard band lower-limit points 
result from comparison with the known limiting results, including those given in the 
above references. 

For M E u = our correlation-function expressions reduce to the well known 
conformal-field general expressions However, for M E u > our expressions provide 
new useful information about the finite-energy spectral properties of the ID Hubbard 
model. 



4.2. INTRODUCTION TO THE FEW-ELECTRON FINITE-ENERGY PROBLEM 

All our applications refer to positive excitation energy u. In this case the general 
correlation function x-si^, can be written in terms of a Lehmann representation as 
follows, 

y oj-ujj,o-ip 

where f3 — > and the operator 0$(k) is the Fourier transform of the general one-electron 
or two-electron physical field <j)$(x, 0) of Eq. (jB.3|) of Appendix B. The j summations 
run now over all available final excited energy eigenstates and Uj t o = [Ei — Egs] are the 
Hubbard-model excitation energies relative to the initial ground state. 

In our applications of the general spectral-function expressions we consider the 
following operators 0#(k): The spin- up one-electron addition operator, the charge 
operator, and the singlet-superconductivity Cooper-pair addition operator, 

&p(fc)=4,t; ^(*) = EE A',a C k+k',a\ 

k' <J=U 

C*ss(^) = X! C fc',| C l-fc',t • ( 22 ) 
k' 

As above, here we use the notations i? = lp for one-electron, d = p for charge, and 
$ = ss for s-wave singlet superconductivity. We are particularly interested in the 
spectral functions associated with the imaginary part of the correlation function (l2T|) , 



lm X 4k, uj) = nJ2\(j\ <%(*0I GS)\ 2 5(lu - u jt0 ) ; = lp, p, ss . (23) 

3 

We consider the case of ground states of spin density m — > 0. In this case the energy 
ujhs — M Eu given in Eq. defines the lower limit of the M th upper Hubbard band. 
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Often the momentum values k l M given in Eq. (jlUJl and associated with the M th upper 
Hubbard band are such that k l M — [tt — I k ] where k is a spectral-function dependent 
momentum. In our applications we consider Hubbard bands generated by dominant 
processes only. According to the results of Ref. [T^], for the spin-up one-electron 
addition operator 0\ p (k) = and charge operator O p (k) = J2k'J2a=i,t c \>,a c k+k',a 
(and singlet-superconductivity operator O ss (k) = J2k' c k' \ c l~-k',i) more than 99% of 
the spectral weight corresponds to excited states such that M Cj _!/ 2 = M = 0, 1 (and 
M C) _ 1/2 = M = 0, 1, 2). States with M Cj _ 1/2 = M > 1 (and M Cj _ 1/2 = M > 2) 
lead to nearly no spectral weight and are omitted below. For electron addition and the 
dynamical structure factor we study the weight distribution for domains of the (k, uj)- 
plane corresponding to small values of (k — k[) and (u — E u ) just above the lower limit 
of the first upper Hubbard band. For the regions corresponding to other momentum 
values and small energy (lo — E u ) these spectral functions vanish because there is no 
spectral weight. In the case of the singlet Cooper pair spectral function we study the 
weight distribution both for domains of the (k, u;)-plane corresponding to small values 
of (k — k[) and (uj — E u ) just above the lower limit of the first upper Hubbard band and 
(k — k l 2 ) and (uj — 2E U ) just above the lower limit of the second upper Hubbard band. 

^From analysis of the energy spectrum of the excited states generated by dominant 
processes, one finds that the M th Hubbard band finite-spectral-weight region of the 
(k, u;)-plane whose lower limit is located at the line uj = M E u is for all values of k 
separated from the next lower-limit line uj = [M + 1] E u of the [M + l] th Hubbard 
band by a region with nearly no spectral weight. The latter region of the (k, o;)-plane 
is out of the range of the excitations generated by dominant processes. This holds 
for all few-electron spectral functions and is illustrated in Fig. 1 for the specific case 
of the one-electron spectral function. For on-site repulsion U = 4t, electronic density 
n = 1/2, and spin density m = the regions of the (k, a;)-plane whose spectral weight 
is generated by dominant processes are shown in the figure by shaded areas. The 
dominant processes also include pseudoparticle particle-hole processes which lead to 
spectral weight both inside and outside but in the close vicinity of the shaded domains 
of Fig. 1. In that figure an extended momentum scheme centered at momentum k = 
(and k = tt) is used for electron removal and the M = lower Hubbard band (and 
M = 1 first upper Hubbard band) . The one-electron addition upper and lower Hubbard 
band spectral functions correspond to values of excitation energy such that uj > E u and 
uj > 0, respectively. Note that for all values of k there is indeed a region with nearly no 
spectral weight located between the finite-weight regions corresponding to the M = 
and M = 1 Hubbard bands. In the figure the one-electron removal spectral function 
corresponds to uj < 0. The figure shows the locations (k = ti — I k F , uj = E u ) for the 
four I = ±1,±3 upper Hubbard band lower- limit points whose weight distribution is 
studied below. There are also two upper Hubbard band lower-limit points for I = ±5 
which are not shown in the figure. 

As in the case of the one-electron spectral weight represented in the figure, it occurs 
for all few-electron spectral functions that in the straight and horizontal lower limit 
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Figure 1. Regions of the (k, w)-plane where accordingly to the restrictions associated 
with the dominant holon and spinon microscopic processes studied in reference 
there is a significant amount of one-electron spectral weight for on-site repulsion 
U = At, electronic density n = 1/2, and spin density m = 0. Note that an extended 
momentum scheme is used. The one-electron addition upper and lower Hubbard band 
spectral functions correspond to values of excitation energy such that lo > E u and 
< lo < E u , respectively. 

upper Hubbard band lines (k, uo = M E u ) of the (k, c<j)-plane there is finite spectral 
weight for a set of discrete momentum values only, as mentioned above. These discrete 
momentum values coincide with the momenta k l M given in Eq. (JTDJ). The four momenta 
k = 7r — / kF where / = ±1, ±3 that are shown in Fig. 1 are examples of such momentum 
values. Each of these momentum values correspond to a different reduced J-CPHS 
ensemble subspace of the CPHS ensemble space spanned by states of rotated electron 
double occupation M and with an occupancy of M — M Cj _i/ 2 —1/2 holons. Our scheme 
provides the weight distribution for small energies (u — M E u ) corresponding to regions 
of the (k, cj)-plane just above such straight and horizontal lower-limit upper Hubbard 
band lines. For the straight and horizontal lines parallel to the lower-limit lines and 
corresponding to values of u such that {u — M E u ) is small, there is finite spectral weight 
for small momentum domains centered around the set of momentum values k = k l M . 

The correlation-function expressions obtained from standard conformal-field theory 
|19[l2~n] correspond to a specific limit of our scheme such that M = 0. The corresponding 
low-energy edges result from ground-state transitions to M C; _y 2 = M = excited states. 
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In figure 1 this corresponds to the finite-weight regions in the vicinity of the three points 
(k = k F , lu — 0), {k = 3k F , to — 0), and {k = 5k F , lu — 0). In the following applications 
of our finite-energy method we do not consider weight distributions resulting from such 
low-energy transitions, which can be evaluated by means of two-component conformal- 
field theory [2H] . Our studies are restricted to the weight distribution of each of the above 
spectral functions in the vicinity of finite-energy points associated creation of —1/2 Yang 
holons and/or cv pseudoparticles of bare momentum q « = ±[71"— 2k F \. The weight 
distributions resulting from creation of cv pseudoparticles at other bare-momentum 
values appear at finite energy uj such that (uj — uj hs ) is finite. These weight distributions 
cannot be evaluated by our method and are studied elsewhere (01 EH EE] . Below we find 
that several critical exponents which control the weight distribution in the vicinity of the 
first and second upper Hubbard band lower-limit points equal corresponding exponents 
that control the low-energy weight distributions studied in Ref. |2Hj- However, we note 
that in the figures that reference these exponents were plotted for finite values of the 
magnetic field whereas our study refers to zero field and magnetization. 

Transitions from the ground state to excited states belonging to the same J-CPHS 
ensemble subspace are of the same type, whereas transitions of different type refer to 
excited states belonging to different J-CPHS ensemble subspaces. In some cases we 
extend this concept to pairs of J-CPHS ensemble subspaces with the same absolute 
momentum value for k l M . 

4.3. ONE-ELECTRON ADDITION UPPER-HUBBARD BAND WEIGHT 
DISTRIBUTION 

We start by considering addition of a spin-up electron. According to Eq. (fT2"|) and 
Eqs. (51) and (52) of Ref. [7 with AL Sj _ 1/2 = and AN S „ = for v > 1, in the 
case of the M c< _i/ 2 = M = 1 type of transitions one has that AN C = AN S = —1 
and AJ S = ±1/2. There are two one — 1/2 Yang holon types of transitions and six 
one-g = ±g° 1 = ±[n — 2k F ] cl pseudoparticle types of transitions. For m — > the 
number of momentum values k = k[ associated with dominant processes such that 
uj = E u is in the present case six. These momentum values read n — I k F where 
/ = ±1, ±3, ±5. Indeed, two one — 1/2 Yang holon type of transitions and two of 
the six one-g = ±g° 1 = ±[n — 2k F ] cl pseudoparticle type of transitions have the same 
momentum and energy spectrum. Thus, they contribute to the same upper Hubbard 
band lower-limit weight distribution located around the points (k = tt k F , uj = E u ). 
These are the two one — 1/2 Yang holon types of transitions such that L C) _i/2 = 1, 
AJ C = 0, and AJ S = ±1/2 and the two one-g = ±g° !/ = ±[ir — 2k F ] cl pseudoparticle 
types of transitions such that iV cl = 1, AJ C = ±1/2, J cl = ±1/2, and AJ S = ±1/2. 
Two other upper Hubbard band lower-limit points are (k = tc =1= 3k F ,uj = E u ). The 
corresponding spectral weight is generated by types of transitions such that iV cl = 1, 
AJ C = ±1/2, J c i = ±1/2, and AJ S = =Fl/2. Below we also study the weight 
distribution in the vicinity of these points. In their vicinity there is a smaller amount of 
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spectral weight than in the vicinity of the above two former points. Finally, the weight 
located in the vicinity of the remaining two points is even smaller. It is generated by 
transitions such that N c i = 1, AJ C = ±1/2, J c \ = =Fl/2, and AJ S = ±1/2. These two 
points are located at (k = tt ^= 5k F ,u = E u ). In the following we study the weight 
distribution around the above six upper Hubbard band lower-limit points located at 
(k — 7T — I k F , lu = E u ) where I = ±1, ±3, ±5. 

Let us consider the limit m — > where the spin-up and spin-down one-electron 
spectral functions have the same form. Use of the general expression ()17j) leads 
for excitation energy u just above E u and momentum values k = it — I kp where 
/ = ±1, ±3, ±5 to, 



To reach this result we first evaluated the expression for general values of m in terms 
of the parameters (jA.2j) of Appendix A. Use of the limiting values for these parameters 
given in the same Appendix then leads to expression (|25j). The same procedure is 
followed in the calculation of the m — > weight distribution of the dynamical structure 
factor and singlet Cooper pair spectral function given below. Note that in the present 
case all the transitions leading to the weight distributions ()24j) involve creation of both 
a c pseudoparticle hole and a s pseudoparticle hole. Note also that the dependence on 
n and U / 1 of the exponent (|25|) occurs through the parameter K p only, as mentioned in 
the previous section. For m = this is a general property of the critical exponents that 
control the weight distribution in the vicinity of the upper-Hubbard bands lower limit. 
In contrast, the expressions of the exponents that control the weight distribution in the 
vicinity of the branch lines studied in Ref. JU] involve the momentum- dependent phase 
shifts defined in Ref. jHJ. 

The exponent (f23J) is a function of U/t which for electronic densities n such that 
< n < 1 changes from Q Z ' — ► {I 2 — 5)/4 = —1, 1, 5 for I = ±1, ±3, ±5, respectively, 
as U/t -> to ('J -> (I 2 - 8)/8 = -7/8, 1/8, 17/8 for / = ±1, ±3, ±5, respectively, as 
U/t — > oo. Note that as the critical exponent f!25l) approaches —1 for I = ±1 as U/t — > 
the spectral function behaves as given in Eq. (|20|). The exponents ()25|) are plotted in 
Figs. 2 (a) and (b) for k = n =p kp and k = tt =p 3kp, respectively, as a function of the 
electronic density n and for different values of U/t. The I = ±5 exponent is not plotted 
in the figure. This exponent is larger than two and corresponds to regions of very little 
spectral weight. The exponents plotted in Fig. 2 are monotonous functions of n. The 
exponent (l p is always negative and such that —1 < (l p < —7/8 and is associated with a 
weight-distribution singularity. It increases for increasing values of U/t. For finite values 
of U / 1 it is a function of the electronic density n with a minimum for an intermediate 
value of n. The exponent (f is always positive and such that 1/8 < < 1 and is 



ah 

Imxip(vr -lk F , u)oc(u- E u )^p ; 
In this case the general exponent (|T9|) is given by, 



/ = ±1, ±3, ±5. 



(24) 




I = ±1, ±3, ±5. 



(25) 
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associated with a weight-distribution edge. The exponent decreases for increasing 
values of U/t. For finite values of U/t it is a function of the electronic density n with a 
maximum for an intermediate value of n. 

For all values of U/t there is more spectral weight around the lower-limit upper 
Hubbard band points located at {k = tt =f kp, uj = E u ) than around the points located 
at (k = tt =1= 3kp, <jj = E u ). Moreover, the amount of weight in the vicinity of the points 
located at (k = n^f3kp, uj = E u ) is larger than that of the weight located in the vicinity 
of the points (k = tt 5kp, to = E u ). For U/t — * the spectral weight located in the 
vicinity of the points {k = tt =F 3k F , u = E u ) and {k = ir =F 5k F , uj = E u ) disappears. 
As mentioned above, the value of the critical exponent plotted in Fig. 2 (a) is —1 for 
U/t —>■ 0. Thus, in this limit the spectral function behaves as given in expression Eq. 
(I20|) . Importantly, in that limit the energy E u of the upper Hubbard band lower limit 
equals the non-interacting electronic spectrum at (k = tt =p kp, uj = E u ). It follows that 
our general U/t weight distribution leads to the correct non- interacting spectrum in the 
U/t — > limit. In this paper we study the weight distribution in the vicinity of points 
located in the lower limit of the Hubbard bands. The branch-line studies of Ref. ^U] 
reveal that the lower-Hubbard band and upper-Hubbard band spectral weights become 
connected at k = tt — kp as U/t — > 0. For finite values of U/t, the finite-weight regions 
of these two bands are separated by a region out of the range of dominant processes. 
This is confirmed by analysis of Fig. 1. 

Our above finite-energy expressions were derived for the metallic phase 
corresponding to electronic densities < n < 1. In the limit n — > 1 and finite values 
of U/t expressions (J2"5j) lead to the correct n = 1 results. In that limit and for m — > 
the upper Hubbard band lower-limit points (k — tt ^ kp, uj = E u ) coincide with the 
lower-limit points (k = tt 3kp, uj = E u ), and (k = 7i =p 5kp, uj = E u ). Thus, in this 
case the smallest of the three exponents (J23j) . such that / = 1, controls the leading-order 
weight-distribution term 

In the limit n — > 1 the equality of the exponent that controls the electron removal 
lower-Hubbard band and the I = ±1 finite-energy electron addition exponent (|2"5|) is 
required by particle-hole symmetry. Our results reveal that in the metallic phase the 
I = ±1, I — ±3, and / = ±5 finite-energy exponents (}25|) also equal the exponents that 
control the weight distribution in the vicinity of the three zero-energy points located at 
[k = ±k F , uj = 0), {k = ±3k F , uj = 0), and {k = ±5k F , uj = 0), respectively. (The 
points (k = kp, uj = 0), (k = 3kp, uj = 0), and (k = 5kp, uj = 0) are shown in Fig. 1.) 
For / = ±1 the exponent expression (|23j) was already found in Ref. |20j . It corresponds 
to the H — > expressions of the one spin-up and spin-down electron exponents given in 
Table I of that reference, where H is the magnetic field. 

There are not many previous results about the excitation energy and momentum 
dependence of the weight distribution in the vicinity of the UHB lower-limit points 
considered above. Conformal-field theory does not apply to finite energy. The method 
used in Ref. jEj refers to U/t — > oo where the UHB lower limit corresponds to 
uj = E u = oo. Thus, the UHB was not considered in the finite-energy studies of that 



One and Two- 



20 





0.2 0.4 0.6 0.8 1 

n 

Figure 2. The upper Hubbard band one-electron exponents (a) £i p and (b) £f p given 
in Eq. (|25(l as a function of the electronic density n for U = 1.5t, J7 = 4t, [/ = lOi, 
and t7 — > oo (horizontal lines). 
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reference. The point (k = n — kp, uj = E u ) corresponds to (k = ir/2, (uj — jj)/t ~ 0.5) 
and (k = 7vr/12 « 0.58 tt, (u-fi)/t w 1) in Figs. 3 and 4 of Ref . [T3], respectively. The 
spectral weight of Fig. 3 is for half filling and U/t = 4. (The figure zero-energy is the 
middle of the Mott-Hubbard gap.) The weight of Fig. 4 is for electronic density n = 5/6 
and U/t = 4. These weight distributions were obtained by numerical calculations 
based on a combination of exact diagonalizations of finite clusters with strong-coupling 
perturbation theory and are consistent with our results. Unfortunately, a quantitative 
comparison is not possible because the method used in Ref. (El does not provide 
accurate information about the weight-distribution dependence on the excitation energy 
and momentum. 

44. THE DYNAMICAL STRUCTURE FACTOR UPPER-HUBBARD BAND 
WEIGHT DISTRIBUTION 

According to Eq. JX2J) and Eqs. (51) and (52) of Ref. |ZJ with AL s ^ 1/2 = and 
AN sl/ = for v > 1, for M c _!/ 2 = M = 1 types of transitions one has in this case that 
AN C = —2, AN S = —1, and AJ S = 0. These transitions generate spectral weight in the 
vicinity of two upper Hubbard band lower-limit points. There are two one — 1/2 Yang 
holon types of transitions and two one-g = ±g° x = ±[n — 2kp] cl pseudoparticle types 
of transitions which contribute to these weight distributions. The former two types of 
transitions correspond to L Ci _i/ 2 = 1 and AJ C = ±1/2. The latter two types correspond 
to N c i = 1, AJ C = 0, and J c \ = =Fl/2. The weight distribution in the vicinity of the 
two upper Hubbard band lower-limit points located at (k = it =p 2kp, uj = E u ) results 
from transitions belonging to these two types. Use of the general expression (fTTjl leads 
in the limit m — ► to the following expression for the dynamical structure factor for 
excitation energy uj such that (u — E u ) is small and momentum values k = ir =p 2k p , 



^From use of Eq. ()19|). we find the following expression for the critical exponent ( p in 
the limit m — > 0, 



The exponent ()27)1 is plotted in Fig. 3 as a function of the electronic density n for 
different values of U/t. This exponent is a function of U/t which for electronic densities 
densities n such that < n < 1 and spin density m = changes from Q p — > as 
U jt —>■ to (, p — > 1/2 as U/t — > oo. For the n = 1 Mott-Hubbard insulator phase the 
exponent reads 1/2 for all finite values of U/t. The exponent ( p is always positive and 
such that < ( p < 1/2 and is associated with a weight-distribution edge. It increases 
for increasing values of U/t. For finite values of U/t it is a function of the electronic 
density n with a minimum for an intermediate value of n. 

In contrast to the above one-electron problem, the spectral weight associated with 
the dynamical structure factor upper Hubbard band vanishes in the limit U/t — ► oo. For 
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(26) 




(27) 
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Figure 3. The exponents £ p = Q s Q27JI that control the weight distribution of the 
upper Hubbard band dynamical structure factor and second upper Hubbard band 
singlet Cooper pair spectral function, respectively, as a function of the electronic 
density n for U = 1M, U = it, U = lOi, and U — >• oo. 



the n = 1 Mott-Hubbard transition, the whole dynamical structure factor vanishes as 
U/t —>■ oo. This behavior is related to the vanishing of the kinetic energy as U/t — ► oo 
|21j . The dynamical structure factor upper Hubbard band exponent (|27jl had not been 
studied until now. It controls an onset of spectral weight whose u derivative is infinite 
at the lower limit of the upper Hubbard band uj = E u and raises upwards according to 
the power law (|26|) . 

In general the momentum values k = 7c^f2kp are finite and thus in the metallic phase 
the weight distribution (|2*Bj) is not related to the zero-momentum frequency dependent 
optical conductivity. However, in the limit n — > 1 these momentum values vanish. 
Thus, for the particular case of the Mott-Hubbard insulator one can use the following 
relation between lmx P (k, uj) and the regular part of the frequency dependent optical 
conductivity a re9 (uj), 

Re^o;)oclim ^# " } , (28) 
fc-^o k 2 

to find the Recr re9 (u;) weight distribution in the vicinity of the upper Hubbard band 
lower-limit point located at (k = 0, uj = Eu = Emh)- Use of this relation reveals 
that the exponent ()27j) controls the frequency dependence of that weight distribution. 
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Combination of Eqs. (|26|) and (|28|) for n — > 1 and E. t 



u 



E MH leads to, 



Rea reg (u) oc (u - E MH ) 



1/2 



(29) 



This expression applies to all finite values of U/t. 

In contrast to the upper Hubbard band dynamical structure factor expression (j2Bj) 
which refers to k = n =f 2k F and values of lj such that (u — E u ) is small, the frequency 
dependence (J2HJ) of the k = optical conductivity in the vicinity of the upper Hubbard 
band lower-limit point was studied by other methods ^J|^j. In reference ^3] that zero- 
momentum conductivity edge was studied for the metallic phase, whose critical exponent 
is different from the exponent ()27|) . However, in the limit n — > 1 the zero-momentum 
density dependent exponent found in such a reference provides the correct half-filling 
exponent. In spite of the different method used in the evaluation of that exponent, the 
half-filling expression derived in Ref. coincides with that of Eq. (J2U|) . Moreover, the 
same optical conductivity edge was studied in Ref. [THj for both small and large values 
of U/t. These investigations also lead to the same critical exponent 1/2 for both these 
limits, in agreement with expression (}29|) which is valid for all finite values of U/t. As 
the conductivity exponent —1/2 is characteristic of a finite-energy semiconductor edge, 
also the conductivity exponent 1/2 is now believed to be an universal signature of the 
Mott-Hubbard insulator [TIHT5] . 

4.5. SPIN-SINGLET COOPER-PAIR FIRST AND SECOND UPPER-HUBBARD 
BAND WEIGHT DISTRIBUTIONS 

The energy E u becomes small for low values of U/t and electronic densities n in the 
vicinity of one. Thus, it is interesting to clarify whether there are singular spectral 
features in the Cooper-pair spectral function at the upper-Hubbard bands lower limit. 
Indeed, for the above values of U/t and electronic density such singular features could 
lead to a superconductivity instability for a system of weakly coupled Hubbard chains. 
Unfortunately, the only singular feature found below is a 5 function peak located at 
the point (k = tt, uj = E u ). That isolated peak results from the ^-pairing mechanism 
[T7] . As the electronic density n approaches one the weight of that peak vanishes as 
(1 — n). While this spectral structure cannot lead to a superconductivity instability for 
the coupled-chain system, other singular features for excitation energy above the first 
upper Hubbard band lower limit could exist, as further discussed in Sec. V. 

According to Eqs. JEJ) and Eqs. (51) and (52) of Ref. with AL s ^ 1/2 = and 
AN SU = for v > 1, in the case of the M c ^x/2 = M = 1 type of transitions there 
is for the spin-singlet Cooper-pair spectral function a one — 1/2 Yang holon type of 
transition and three one-g = ±g° x = ±[n — 2kp] cl pseudoparticle types of transitions. 
However, we find that the one-g = ±g° x = ±[n — 2k F ] cl pseudoparticle transition 
whose excitation momentum and energy is the same as the momentum and energy of 
the one — 1/2 Yang holon type of transition does not contribute to the singlet Cooper 
pair spectral function. Thus, only the one — 1/2 Yang holon type of transition and two 
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of the three one-g = ±g° x = ±[n — 2k F ] cl pseudoparticle types of transitions contribute 
to weight-distribution features. We find that the weight-distribution feature generated 
by the one — 1/2 Yang holon transition is a single S peak located in the first upper 
Hubbard band lower limit at k = n and uo = E u . The matrix element between the 
corresponding cl pseudoparticle excited state with the same momentum k = tt and 
excitation energy u = E u and the ground state vanishes in this case. Therefore, this 
singlet Cooper pair spectral function feature results from the one — 1/2 Yang holon 
excited state only. There are also two first upper Hubbard band lower-limit points 
located at {k = TT^f4kp, uo = E u ) which are generated by one-g = ±<7°i = ^t 71 " — 2k f] cl 
pseudoparticle transitions. Below we also study the weight distribution in the vicinity 
of these first Hubbard band lower-limit points. 

For the one — 1/2 Yang holon transition we find that AN C = AN S = AJ C = AJ S = 
0. Thus, this is an example of a transition where the quantity ()19j) reads ( ss = —2 
and is not a critical exponent. Instead, the correlation function at that point of the 
(k, u;)-plane is of the form given in Eq. (|2T)|) for all values of U/t and reads, 

C 

Xss{rr, u) = =r — ~3 > ( 30 ) 

uj - E u - ifj 

where (3 — > and C p > is a constant. It follows that the associated singlet Cooper 

pair spectral-function weight distribution is indeed of the type given in Eq. (|2T)|) and 

reads, 

Imx ss (vr,u;) = rrC p 5(uj - E u ) . (31) 

The weight feature (|31|) can be obtained by direct evaluation of the matrix elements 
on the right-hand side of expression (}2Tj) for $ = ss. That procedure confirms in this 
particular case the validity of the results obtained by our method. Moreover, direct 
calculation of the matrix elements reveals that for k = tt the expressions ()30|) and (|3H) 
do not refer to values of u such that (u — E u ) is small only, but to all values of excitation 
energy u. Direct evaluation of the matrix elements also reveals that the constant C p of 
Eqs. (|3Uj) and (}3*Tj) reads C p = [N a — N] and vanishes for the Mott-Hubbard insulator. 

In order to evaluate the matrix elements of expression fl21(l for $ = ss and in the 
particular case of the momentum value k = tt, we note that the off diagonal generator 
SI of the ?7-spin SU(2) algebra given in Eq. (9) of Ref. |Hj can be rewritten as follows, 

= X! c k',i c I-fe',T • (32) 
k' 

Application of that generator onto a ground state produces an excited state with a —1/2 
Yang holon. The generator of the above transition corresponds to creation of that —1/2 
Yang holon and involves the initial ground state and this single final excited state only. 
After normalization this excited state can be written as follows, 



, _ J2k' c k',i c t-i '. i 
" V2 " ^ " VK-N 



1) = ^ k *- W A | GS) . (33) 
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The main point is that at k = it the singlet superconductivity operator given in 
Eq. fl21 reads, 

a s (7r) = ^ = £4',|cU, r (34) 
k' 

Use of Eq. ()34|) in Eq. (|21|) reveals that for k = it the correlation function expression 
results from the overlap of the ground state with the excited state (f3~3J) only. Thus, 
for k = 7r only one matrix element between the ground state and the available excited 
states of momentum k = it is finite. That matrix element corresponds to the state 
(jnnj)- It follows that for k = it the singlet superconductivity correlation function and 
corresponding spectral function are indeed given by expressions (JHUj) and (|HTJ) with 
C p = [N a — N] for all values of uj. 

The one — 1/2 Yang holon state (j33|) is associated with the ^-pairing mechanism 
[T7] . That mechanism leads in the metallic phase to the 5-function peak (|3~Tj) in the 
singlet Cooper pair spectral function. Such a state possesses off-diagonal long-range 
order |17j . The corresponding 5 peak is located at finite excitation energy. The 5 peak 
(f3Tj) is absent in the case of the Mott-Hubbard insulator where C p = [N a — N] — > 0. This 
vanishing is required by a half-filling selection rule, since the ground-state 7] spin value 
of the Mott-Hubbard insulator reads S c = 0. Note that the present transition is a mere 
rotation in 7/-spin space which conserves S c and leads to AS C = and AS* = +1. Thus, 
only for initial ground states such that S c > 1 can this transition occur. This excludes 
both the Mott-Hubbard insulator ground state such that S c = and the one-hole doped 
Mott-Hubbard insulator ground state such that 5* c = 1/2. 

Next, we consider the weight distribution in the vicinity of the first upper Hubbard 
band lower-limit points generated by other M C) _i/2 = M = 1 types of transitions. These 
transitions correspond to a region of little spectral weight and are such that iV c i = 1, 
AJ C = ±1/2, J cl = =Fl/2, and AiV c = AN S = AJ S = 0. Use of the general expression 
fTTj) leads in the limit m — >• to the following expression for the spin singlet Cooper 
pair spectral function for momentum values k = it =p 4k p and excitation energy uo such 
that {uo — E u ) is small, 

Imxs s (vr =F Ak F , u) oc (a; - E^ . (35) 

This weight-distribution edge does not occur for the Mott-Hubbard insulator. From use 
of Eq. (fTT?|) we find the following expression for the critical exponent ( ss in the limit 
m — > 0, 



C SS = -2[1-[^2K P Y\. (36) 

This exponent is a function of U / 1 which for all electronic densities changes from ( ss — > 2 
as U/t — > to ( ss — >• as U/t — >• oo. Such a dependence on U/t is continuous. (Since 
there is not much spectral weight in the vicinity of this edge feature, we do not plot the 
positive exponent (jHEJ)-) 

We close our study by considering the weight distribution in the vicinity of the 
M c -i/2 = M = 2 second upper Hubbard band lower-limit points located at excitation 



One and Two-Electron Spectral Function Expressions. 



26 



energy u = 2E U . We note that for the Mott-Hubbard insulator there is no first upper 
Hubbard band. In this case there is an energy gap which equals twice the Mott-Hubbard 
gap and separates the singlet Cooper pair removal and addition spectral functions. The 
latter function corresponds to the present second upper Hubbard band in the limit 



The M Cj _i/2 = M = 2 types of transitions are such that AiV c = —2, AN S — — 1, 
and AJ S = 0. Interestingly, there are five different types of transitions which contribute 
to the same weight distribution in the vicinity of two points and lead to the same critical 
exponent. These five types of transitions are such that (a) L c _i/ 2 = 2 and AJ C = ±1/2, 



(b) Lc-i/a = 1, N cl = 1, AJ C = 0, J cl = ±1/2, (c) N cl = 2, AJ C = ±1/2, J cl = 0, 
(d) N cl = 2, AJ C = ±1/2, J cl = ±1, and (e) N c>2 = 1, AJ C = 0, J cl = ±1/2. All 



transitions of these types contribute to the weight distribution around the two second 
upper Hubbard band lower-limit points located at (k = ±2kp, u> = 2E U ). In addition, 
there is a sixth type of transition which leads to two second upper Hubbard band lower- 
limit points located at (k = ±6kp, oo = 2E U ). This type of transition is such that 
N c i = 2, AJ C = ±1/2, and J cl = =Fl/2. There is nearly no spectral weight in the 
vicinity of these two latter points. 

Use of the general expression (JT7j) leads in the limit m — > to the following 
expression for the spin singlet Cooper pair spectral function for excitation energy uj 
such that (u — 2E U ) is small and momentum values k — k\ — 1 2kp where I = ±1, ±3, 



^From use of Eq. (|19|). we find the following expression for the critical exponent Ci'J m 
the limit m — > 0, 



Note that for / = ±1 this exponent equals the dynamical structure factor exponent 
( p given in Eq. (|77j) . Thus, it changes from (l s — > as U/t — > to (l s — > 1/2 as 
U/t — > oo. It is plotted in Fig. 3 as a function of the electronic density n for different 
values of U/t. On the other hand, the exponent is much larger and changes from 
CL ~* ^ as U/t — > to (l s — > 9/2 as U/t — > oo. This exponent corresponds to a 
region of very little spectral weight. In the limit n — > 1 the four second upper Hubbard 
lower-limit points located at {k = ±2kp, to = 2E U ) and (k = ±6kp, uj = 2E U ) become 
the same single point. Thus, in this case the smaller exponent C,l s corresponds to the 
dominant contribution and controls the weight distribution in the vicinity of the lower- 
limit point located at {k = ir, oj = 2Eu = 2Emh)- For the n = 1 Mott-Hubbard 
insulator phase the exponent (l s reads 1/2 for all finite values of U/t. It equals the 
exponent which controls the weight distribution around the zero-energy point located 
at (k — 7T, u — 0). Interestingly, comparison with the results of Ref. J2D] reveals 
that also in the metallic phase the exponents ()38j) for I = ±1 and I = ±3 equal the 
corresponding exponent that control the low-energy weight distribution of the singlet 
Cooper pair removal spectral function in the vicinity of the zero-energy points located 



lm.Xs S (l2k F , lo) oc (u - 2E U ) 



±1, ±3. 



(37) 




±1, ±3. 



(38) 
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at (k = ±2kp, uo = 0) and (k = ±6kp, u — 0), respectively. However, in the metallic 
phase there is more spectral weight in the vicinity of the latter zero-energy points than 
in the vicinity of the corresponding second upper Hubbard band lower-limit points. 

Finally, we note that if one extends the present study to finite-energy excited states 
with finite occupancies of —1/2 HL spinons and/or sv pseudoparticles belonging to 
v > 1 branches, similar features are obtained for gapped spin excitations. For instance, 
in the case of the spin-down triplet superconductivity spectral function a 5-function 
peak similar to (pTTj) is obtained at k = ir and u = E u + fi H. This peak is generated by 
a transition from the ground state to a single finite-energy excited state with one —1/2 
Yang ho Ion and one —1/2 HL spinon. 

5. CONCLUDING REMARKS 

In this paper we derived general few-electron spectral function expressions for the ID 
Hubbard model in the vicinity of the lower limit of the upper Hubbard bands. For the 
one electron addition and dynamical structure factor we studied the weight distribution 
in the vicinity of the lower limit of the first upper Hubbard band. In the case of the 
singlet Cooper pair spectral function we considered the same problem in the vicinity 
of the lower limits of both the first and second upper bands. The weight of these 
Hubbard bands is generated by dominant holon and spinon processes which amount to 
more than 99% of the spectral weight [Ej. Our general expressions were obtained by 
combination of symmetries of the ID Hubbard model, such as the holon and spinon 
numbers conservation laws, with the finite-energy holon and spinon description of the 
quantum problem recently introduced in Refs. [HI IE] • 

The results obtained in this paper provide physically interesting and useful 
information about the finite-energy spectral properties of the present many-electron ID 
quantum liquid and are complementary to the branch line studies of Refs. j^J EH 
In the case of the one-electron problem the combination of our results with those of Ref. 
[TU] provides all finite-energy spectral-weight singularities. Since the singular spectral 
features observed in quasi-lD metals agree quantitatively with the model singular branch 
lines, our results are of interest for the further understanding of the unusual spectral 
properties observed in low-dimensional materials [21 IU El El GUI- Moreover, from the 
finite-energy weight-distribution found for the dynamical structure factor, we checked 
that in the limit of half filling our results lead for all finite values of U/t to the 
expected Mott-Hubbard insulator exponent 1/2 for the onset of the finite-frequency 
optical conductivity absorption. Unfortunately, we found no singular features at the 
the upper-Hubbard bands lower limit of the spin-singlet Copper-pair spectral function 
other than the expected 77-pairing ^-function peak. However, it could be that there are 
branch lines just above such a lower limit which show singular spectral features. Since 
the energy width of the first Hubbard band is small for low values of U / 1 and electronic 
densities n in the vicinity of one, such singular features could lead to an instability for a 
system of weakly coupled Hubbard chains. We thus suggest that this problem is further 
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studied by the branch-line method used in Ref. ^U] for the one-electron problem. 
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Appendix A. THE FINITE-SIZE ENERGY SPECTRUM QUANTITIES 

In this Appendix we show that the functional 2 has the form given in Eq. ()15|) 
and discuss issues related to the vanishing of the cv finite-size energy spectrum for 
q = ±ql = ±[n-2k F ]. 

By direct use of the results of Ref. [Hj, we find that for the excited states that span 
the reduced J-CPHS ensemble subspaces H re d the functional (pMj) can be written as, 

*K=[<> E t^ + ^Cy 

a'=c,s u=l 

+ E ^AJ^-E^J^] . (A.l) 

a'=c,s v=\ 

Here the parameters a , and cv can expressed in terms of the two-pseudofermion 
phase shifts defined in Ref. jH] as follows, 

l=±l 

l=±l 

In the limit m — > these parameters read £° c = l/^j2K p , = 0, £° c = —l/y/2, 

e s = y/2, e° = o, e° = o, # c = = e., e = °> £. = ws, 



£c,ck = Y^-^P) an d ^s,o/ = 0- Here K p is the parameter defined in Ref. |22] which is 
such that K p — > 1 as — ► and — > 1/2 as U/t — > oo. 

The creation of —1/2 Yang holons does not lead to any contribution to the value of 
the functional (|A.1|) . However, the creation of cv pseudoparticles leads to contributions 
through the terms iY^=i £.a,cuN C u/2 and — Y^Li ^x CV Jcv of that functional. On the 
other hand, the symmetry discussed in the text below Eq. (JT5|) requires that the values 
of the phase-shift parameters cv of Eq. (|A.2jl are such that the functional (jA.lj) 
simplifies to (fTK|) . Indeed, by manipulation of the integral equations of Ref. [Hj that 
define the two-pseudofermion phase shifts, we find the following relation between the 
parameters cv and ^ c defined in Eq. (IA.2J) . 



€ 3 a,cv = 3&,c> a = c,s; u = 1, 2, 3, ... ; j = 0, 1 . (A.3) 

Use of this relation in Eq. (|A.1J) leads to (fT5J) . 

Moreover, the same symmetry is behind the vanishing of the finite-size energy 
contributions of the = ±[n — 2k f] cv pseudoparticles. Indeed, by careful analysis 
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of the general energy spectrum AEql we find that these contributions vanish provided 
that the quantity, 



2 At 



cv 



'tAq cv + ; u=l,2, ... ; i = ±1 , (A.4) 



-2tt V w A a 

also vanishes. Here Q C v{q) is the functional defined by Eq. (73) of Ref. |Sj. Use of the 
equality Q C u(Lq%)/N a = —iAq® u given in Eq. (126) of the same reference in expression 
(jA.4j) confirms that 2 A L CV vanishes for the excited states that span the reduced subspace. 

These properties apply only to cv pseudoparticles such that q — > ±[K—2kp]. Indeed, 
these quantum objects become non-interacting and localized in that limit. This behavior 
is related to the fact the —1/2 Yang holon is also non-interacting and invariant under 
the electron - rotated-electron unitary transformation However, in general a cv 
pseudoparticle is different from the rotated cv pseudoparticle. The only exception is 
precisely for bare momentum values q such that q — > ±q% = ±[7r — 2kp]. The cv 
pseudoparticles become non-interacting in that limit only. Indeed, the energy associated 
with creation of a —1/2 Yang holon is given by E u . A cv pseudoparticle is a composite 
quantum object of v —1/2 holons and v +1/2 holons. Creation of a +1/2 holon is a 
process which requires no energy. The amount of energy required for creation of a cv 
pseudoparticle is v E u + e° cv (q). Thus, the energy per —1/2 holon is E u + t° cv (q) jv ^ E u . 
However, one has that e° cu {q) — > as q — > ±g°, = ±[tc — 2k F ] and thus in this limit 
the energy per —1/2 holon becomes E u and equals that of a non-interacting —1/2 Yang 
holon. 



Appendix B. THE FINITE-ENERGY CORRELATION FUNCTION 

Here we show that the asymptotic expansion of the correlation function 
(GS\ (f>$(x, t) 0,?(O, 0) | GS) is of the form given in Eq. (flfi|) . The operator (j)^ L (x, t) acts 
and is defined in a reduced J-CPHS ensemble subspace associated with small values of 
the excitation energy AE GL = (uj — u HS ) and excitation momentum AP GL — (k — k l M ). 
In the Heisenberg description and in that reduced Hilbert subspace the space coordinate 
x and time t dependent physical field 0f L (x, t) can be expressed relative to its x = 
and t — expression 0#(O, 0) as, 

j) G d L (x, t) = e -i[-HGL-t-.P G L-x) q) e i{:H GL :t-.P GL :x] _ (g J) 

For the Hamiltonian : Hql '■ of Eq. © the low-energy energy spectrum ()13|) and 
momentum spectrum APql = jf E a =c lS J2 L =±i 1 [ A L a + N*£ L ] are conformal invariant 
and correspond to a two-component conformal-field theory associated with the c and 
s pseudoparticle excitation branches 120] • The effect of the cv pseudoparticles is 
merely to shift the c pseudoparticle current number deviations AJ C by — Y^=\ Jcv 
Fortunately, such an effect does not affect the conformal invariance of the low-energy 
energy spectrum and momentum spectrum. Thus, the asymptotic expression for the 
low-energy correlation function of the physical field 0^ L (x, t) is for electronic and spin 
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densities such that < n < 1 and < m < n, respectively, of the following general form 

[111201, 

(GS\jg L (x,t)$$ L (p,0)\GS)<x J] II 7Z _ * f ^ • ( B -2) 

a=c, s L =±l \ X LVa l ) a 

We recall that the ID Hubbard model : H : and Hamiltonian : Hql '■ of Eq. (JHJ) 
refer to the same ground state. Thus, the field <f>#(x, t) is such that, 

Mx, t) = < ■ " ' :/ ' : r MO, 0) < > Jhi ip:x , (B.3) 

where : P : is the momentum operator of Eq. (JHJ). Since the Hamiltonians : H : and 
: Hql '■ °f Eq. (JHJ) have the same energy eigenstates, there is a one-to-one correspondence 
between the correlation function terms (jB.2|) of the Hamiltonian : Hql '■ and those of 
the ID Hubbard model : H :. Transitions to different reduced subspaces 'H re d lead 
to different correlation-function terms. Moreover, transitions to reduced subspaces with 
different values for M E u and k l M lead for the ID Hubbard model to correlation-function 
contributions with different values of the excitation energy uj « M E u and momentum 
k « k l M . Thus, our first choice is the energy value M E u and momentum value k l M 
which our correlation-function asymptotic term refers to. If for the same values of 
M E u and k l M there are several reduced subspaces, we choose that associated with the 
leading-order asymptotic correlation-function term for the Hamiltonian : Hql '■ of Eq. 
(JHJ). That leading-order term is of the form (|B.2J) . Let \ip(j)) with j = 1,2,... denote 
the excited states that span the chosen reduced Hilbert subspace TL re d- Based on the 
commutation relations of the three normal-ordered Hamiltonians related by Eq. (JHJ) and 
three normal-ordered momentum operators involved in Eq. (JHJ) we find that, 



fo( X , t)=e -^HS + -.HcL-.)t-(P + -.PcL-.)x] ^ (0> 0) jl(H H S + -.Hc L -.)t-(Po+-.PcL-.)x] 
— e -i[HHSt-Pox] e -i[:H G L-t-:PGL-x\^ ) jQ Q\^i[Hat t—PGL %] e i[Hnst-Pox] 
= e -^ HS t-P x] ^GL^ t y[H HS t-P x] ^ (B.4) 

where we have used Eqs. (JB.IJ) and ()B.3|) . Equation (JB.4J) confirms that 0^(0, 0) = 
0f L (O, 0). 

Once the ground state has eigenvalue zero both for Hhs and Pq and all excited 
states that span the subspace H re d have for these operators the same eigenvalues M E u 
and k l M , respectively, we find that, 



(GS\fo(x, t) M0, 0)\GS) 

= (GS\e- i[ * Hst - pox] $ L (x,t) e i[flHst - pox] 4>^ L (0, 0)\GS) 

= -£(GS\4>$ L (x, t )^st-P 0X ] |Vi(j)>^(;)l#? L (0, 0)|G5> 
j 

= e *i ME ^-k' M x] (GS\$ L (x, t) 0f L (O, 0)\GS) . (B.5) 



One and Two-Electron Spectral Function Expressions. 
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Here the j summation refers to the above set of excited states which span 

the suitable reduced subspace T~C re d- Finally, the combination of Eqs. (jB.2|) and (jB.5|) 
implies that the leading term in the asymptotic expansion of the correlation function of 
the ID Hubbard model is of the form given in Eq. (|16|) . 
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